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To explore underlying mechanism for the superconducting phase in recent discovered aromatic 
hydrocarbons, we carry out the first-principles calculations on benzene, the basic and the simplest 
unit of the series and examine the structural and phase stability when doped by potassium, K^CeHe 
(x = 1,2,3). We find that K 2 C 6 H 6 with the space group of Pbca is the most stable phase with 
superconducting transition temperature around 6.2 K. Moreover, we argue that all existing hydro¬ 
carbons should have a unified superconducting phase in the same temperature range of 5—7 K, when 
doped by two potassium atoms. Our results indicate that the electron-phonon interaction is enough 
to account for the superconductivity of this unified superconducting phase. 
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Organic superconductors are unique materials with 
a crystal structure made primarily of a complex car¬ 
bon based network. Carbon, an element associated di¬ 
rectly with life, which was postulated to have a high 
transition temperature even above room temperature, 
from a theoretical viewpoint [1]. Organic superconduc¬ 
tors are generally charge-transfer compounds that are 
made up of electron donor molecules such as tetrathia- 
fulvalene (TTF), bis-ethylenedithrio-TTF, (BEDT-TTF, 
abbreviated as ET), and tetramethyltetraselenafulva- 
lene (TMTSF) derivatives, as well as electron acceptor 
molecules such as tetracyanoquinodimethane (TCNQ) or 
fullerides [2-6]. The recent discovery of superconductiv¬ 
ity in polycyclic aromatic hydrocarbons (PAHs) [7-10] 
has led to a great interest in these organic based com¬ 
pounds. Increasing the number of benzene rings, thus 
delocalizing the charge over a larger space, superconduct¬ 
ing transition temperature, T c , could be higher than 30 K 
and thus revitalizing the hope of realizing high-T c super¬ 
conductivity in light element materials. However, there 
are issues to be addresses such as what the supercon¬ 
ducting phase of these hydrocarbons is? how many car¬ 
ries are needed to induce superconductivity? and where 
the doped carriers are distributed in the unit cells, etc. 
Having such information is important not only for un¬ 
derstanding the physical properties and mechanism of 
superconductivity in these materials but also for design¬ 
ing materials with higher T c in technology applications. 

The simplest aromatic hydrocarbon is nothing but the 
benzene, which has six carbons on a ring and six hydro¬ 
gens attached, CeHe- Other PAHs are made of more rings 
of carbons and hydrogens. Being the basic cell of PAH, 
solid benzene stands out as an ideal system to address 
all concerned issues mentioned above. A fundamental 
question is whether such a system could superconduct 
when doped by alkali atoms. If it super conducts, will its 
properties be of common features for all aromatic hydro¬ 


carbons? Thus, explore structures of doped solid benzene 
and examine corresponding phase stability, followed by 
the study of superconductivity, is the task of this work. 

Solid benzene has rich phase diagram with several 
phases as functions of temperature and pressure, whose 
crystallization characteristics has been extensively inves¬ 
tigated [11-17]. Liquid benzene crystallizes at room tem¬ 
perature and about 0.07 GPa, in an orthorhombic phase 
with space group Pbca (phase I), containing four formula 
units (f.u.) per unit cell [11]. At 1.4 GPa, a transition 
occurs from phase I to phase II (P4s2i2), and then at 
4 GPa, phase II transfers to phase III (P2i/c), followed 
by phase IV (11 GPa) and phase V at high pressures. 
Solid benzene is non-metallic and non-superconductive 
in a very large range of pressure. Upon heavy com¬ 
pression, solid benzene was predicted to become a metal 
from insulator in the pressure range of 180 — 200 GPa 
[18], which means that benzene is a possible supercon¬ 
ductor under high pressure. However, extreme pressure 
condition brings the difficulty to experimental fabrica¬ 
tion and observation. Hence, we choose potassium (K) 
to dope solid benzene at ambient pressure to explore 
the realization of metallization and superconductivity. If 
doped solid benzene superconducts at ambient pressure, 
it would be possible to have higher T c at high pressure. 

Although experimental studies on doping K into solid 
benzene have not been performed yet, the interaction be¬ 
tween benzene molecule and K + ion has been observed 
experimentally and studied theoretically by E. Lopez et 
al [19]. Their work implies the feasibility of doping 
potassium into solid benzene. Notably, doping potas¬ 
sium has been shown to be the most effective and com¬ 
mon method to induce superconductivity in PAHs [7-10]: 
all aromatic hydrocarbons were found to become super¬ 
conductor with K dopant. One reason is that K has a 
relatively low melting point as compared to other alkali 
elements, which enables the chemical reaction and crys- 
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FIG. 1: (Color online) Optimized crystal structures of (a) 
solid benzene and (b) K-doped solid benzene, K 2 C 6 H 6 , view¬ 
ing from different directions. Purple balls represent K atoms. 


tallization from liquid states of K and benzene easy to 
occur. 

For doped solid benzene, we firstly perform structure 
optimization to obtain its atomic positions and lattice 
constants. Then we calculate corresponding electronic 
structures and electron-phonon interactions. All struc¬ 
tural optimizations and electronic properties were cal¬ 
culated by using the Vienna ab initio simulation package 
(VASP) [20]. In this work, the projector augmented wave 
(PAW) method [21] was adopted, and the exchange corre¬ 
lation energy was described by the Ceperley-Alder local 
density approximation (LDA) [22] as parameterized by 
Perdew and Zunger [23] , which has been confirmed to be 
adequate by previous studies [24-28]. An energy cutoff 
of 500 eV was used for the plane wave basis sets. The 
Brillouin zone (BZ) was sampled by 4 x 4 x 4 Monkhorst- 
Pack Appoint grids during the optimization, and double 
Appoint in the electronic structural calculations. The con¬ 
vergence thresholds were set as 10 -6 eV in energy and 
0.005 eV/A in force. A conjugate-gradient algorithm was 
used to relax the ions into their instantaneous ground 
state. 

To calculate electron-phonon interaction and relevant 
parameters, we use the QUANTUM-ESPRESSO package 
(QE) [29] with a cutoff energy of 60 and 450 Ry for wave 
functions and charge densities, respectively. Forces and 
stresses for the converged structures were optimized and 
checked to be within the error allowance of the VASP and 
QE codes. Troullier-Martins norm-conserving scheme 
[30] was used to generate the pseudopotentials for C, H, 
and K atoms. 8x8x8 Monkhorst-Pack Appoint grid with 
Gaussian smearing of 0.03 Ry was used for the electron- 
phonon interaction matrix element calculations at 2 x 2 x 2 
g-point mesh. 

As a comparison and a check on accuracy, we start 
our investigation by looking at the solid benzene. Fig¬ 


ure 1(a) shows the optimized crystal structure of solid 
benzene viewing from different directions, with crystal 
lattice constants being a = 7.041 A, b = 8.903 A, and 
c = 6.357 A at zero temperature and the ambient pres¬ 
sure, which are 3-6% less than experimental values at 
78 K [11]. Considering temperature effect, this error be¬ 
tween theoretical prediction and experimental measure¬ 
ment is acceptable. From reported data on PAHs, the 
superconductivity is highly sensitive to both the doping 
concentration and the cation positions in the unit cell. 
To examine both effects, we perform optimization studies 
on doped solid benzene, K^CeHe, by varying potassium 
contents (x = 1, 2, 3) and position of K atoms in the unit 
cell. To check their stability, formation energy, as defined 
by Eformation — Ek x c q Hq E( 3 6 h 6 Ek x , is a good in¬ 
dication. We obtained the results of formation energy as 
0.35, —2.99, and 1.52 eV for x = 1, 2, and 3, respectively. 
Thus, both KiC 6H 6 and K 3 C6H 6 are unstable and the 
K 2 C 6 H 6 is the only stable phase in this series. Its big 
negative formation suggests easier synthesization. 

Figure 1(b) shows the optimized K^CeHg crystal struc¬ 
ture viewing from different directions, with lattice con¬ 
stants being a = 8.605 A, b = 10.705 A, and c = 6.415 A, 
respectively. Comparing with the undoped case, the lat¬ 
tice constants expand largely in the a and b directions. 
The inserting of eight K atoms makes the volume in¬ 
creased by 26.9%. Under K doping, a visible layer shape 
structure is formed, differing from that of the solid ben¬ 
zene. As shown in Fig. 1(b), all benzene rings are ro¬ 
tated and parallel to the froplane under the interaction 
of dopants, and K atoms are intercalated into two lay¬ 
ers of benzene rings, occupying on Wyckoff 8c sites. We 
found that the optimized structure of K 2 C 6 H 6 has tended 
to the higher symmetry of Fmmm from the lower Pbca 
symmetry. Significantly, the layered structure of K 2 C6H 6 
is similar to those of YbCg and CaC6 [31]. This similarity 
indicates the possibility of superconducting, because it is 
known that both YbC6 and CaC6 are superconductors, 
with the T c reaching to 6.5 K and 11.5 K, respectively 
[31, 32]. 

Once the crystal structure is determined, it is ready 
to obtain the electronic band structure. For K 2 C 6 H 6 , we 
show 8 conduction bands (called as CB1, CB2, ..., and 
CB8) along a high symmetry /c-direction, Y — Z — T — 
Y — S — X — U — R, of the Brillium zone in Fig. 2(a). 
Specifically, under either Pbca or Fmmm symmetry, the 
LUMO and LUMO+1 orbitals split into 8 bands from Y 
to Z, then they merge into 4 degenerate bands along the 
Z — T — Y — S — X — U directions, and further into two 
highly degenerated bands along the U—R direction. Elec¬ 
trons transferred from the K atoms to the CeH 6 molecule, 
occupying low energy orbitals and raise the Fermi level. 
Three bands, CB1, CB2, and CB3 are fully occupied by 
six of eight electrons, while the other three, CB4, CB5, 
and CB6 are partially filled by the remaining two elec¬ 
trons with Fermi level passing through. Both CB7 and 























3 



FIG. 2: Electronic band structure of K 2 C6H 6 : (a) energy 
band; (b) Fermi surface; and (c) projected density of states 
(PDOS), where zero energy denotes the Fermi level. 


FIG. 3: Electron-phonon coupling of K 2 C6H 6 : Eliashberg 
spectral function a 2 F(uj) and the electron-phonon coupling 
integral A(a;). 


CB8 are empty. Thus, K 2 C 6 H 6 is metallic and it is in 
the low-spin state. Comparing with the solid benzene, 
the intermolecular distance along benzene rings becomes 
smaller on the 6c-plane of K 2 C 6 H 6 , which enhances the 
intermolecular interaction and broadens the dispersion of 
energy bands. Analyzing the electronic characters near 
the Fermi level, as shown by Fig. 2(b), the hole-like Fermi 
surfaces are formed by CB4 along the Y — Z ^-direction, 
while both CB5 and CB6 contribute the electron-like 
Fermi surfaces around the Y Appoint. Fermi surface nest¬ 
ing is quite clear. Calculation of the projected density 
of states (PDOS) show that C -p y and C-p^ states are 
fully occupied and far below the Fermi level, while C-p^ 
state has a big PDOS around the Fermi level, as shown 
in Fig. 2(c). Most PDOS of K-(s,p x ,p y ,p z ) states are 
above the Fermi level, which indicates that the electrons 
transfer from K to C-p^ orbital. The large ratio of C-p^ 
PDOS over K PDOS shows that C-p^ electronic states 
contribute more to the Fermi surfaces. 

Now we examine whether this simplest aromatic hydro¬ 
carbon could exhibit superconductivity like other PAHs 
wiith long benzene rings [7—10] . Assuming the pair¬ 
ing mechanism is still the electron-phonon interaction 
mediated, we evaluate the Eliashberg spectral function 
a 2 F(uj) and the electron-phonon coupling integral A (a;) 
of K 2 C 6 H 6 as a function of frequency. The results are 
shown in Fig. 3. The vibration in the low frequency range 
of 0 — 250 cm -1 comes from K atoms and intermolecular 
modes, while the vibration with moderate frequencies of 
300 — 1450 cm -1 mainly comes from C atoms. The total 
A is 0.67, which is comparable to that of K 3 C 22 H 14 [25]. 
We also noted that vibrational modes from C atoms con¬ 
tribute 56% to the total A value, and the rest 44% is pro¬ 
vided by the dopant and intermolecular phonon modes. 
Following the modified McMillan equation of Allen and 
Dynes [34], with the obtained values of A (a;) and uoi og 


for K 2 C 6 H 6 , and a typical value of 0.10 for the Coulomb 
pseudopotential p*, we obtain T c = 6.2 K. The choice of 
p* is appropriate for light-element systems [25, 35, 36]. 

We want to point out that the dominant contribution 
to the total A from the vibrational modes of the C atoms 
in K 2 C 6 H 6 shows some similarity with other C6-based 
or low carbon-bearing compounds. In reality, CaC6 was 
discovered to exhibit superconductivity at 11.5 K [31]. 
Superconductivity has also been discovered in interca¬ 
lation compounds of graphite with alkali metals ACg 
(A =K, Rb, or Cs) [37]. The predicted T c of 6.2 K in 
K 2 C6H 6 is comparable to the temperature range of these 
low carbon-based superconductors [31, 37]. 

A comparison of the predicted superconductivity in 
benzene with other aromatic hydrocarbons reveals a com¬ 
mon feature of PAH (Fig. 4). Besides the reported high- 
T c phases with T c = 18 K for K-doped picene ( 1 ^ 022 H 44 ) 
[7], 15 K for K-doped coronene ( K x C 2 /Jii 2 ) [ 8 ], and 33 
K for K-doped l,2:8,9-dibenzpentacene (K^CsoHig) [10], 
multi-superconducting phases in long-benzene-ring hy¬ 
drocarbons were generally observed. There exists a low 
T c phase with almost the same T c between 5 — 7 K in 
these hydrocarbon superconductors discovered so far [7- 
10]: K x Ci 4 Hi 0 has T c 5 K [9]; K x C 22 ^i 4 has T c 7 K 
[7, 8 ]; K^CsoHig has T c of 7.4 K [10], though several 
other transitions with higher T c ’s were observed. More¬ 
over, RbsC 22 Hi 4 and Cai. 5 C 22 Hi 4 both show transition 
temperature around 7 K [7]. T c of 5 — 7 K was re¬ 
ported for K x C 24 Hi 2 [ 8 ], in addition to other three more 
transitions. Superconductivity at 6.4 K was observed in 
Sm 1.25 chrysene [38]. It is important to note that only 
in the range of 5 — 7 K transition temperatures were 
observed in all reported superconducting aromatic hy¬ 
drocarbons with various number of benzene rings, in a 
good agreement with our current prediction for benzene, 
the simplest compound of this series. In fact, the value of 
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FIG. 4: Variation of superconducting transition temperature 
T c with the number of benzene rings for various aromatic 
hydrocarbons. The solid red square is the prediction of this 
work for benzene. Other symbols are taken from experiments 
[7-10, 38]. The dashed line is the guide for eyes. 

the logarithmic average of phonon frequencies uJi og (199.8 
K) for K 2 C 6 H 6 is exactly the same as that of K 3 C 22 H 14 
[25], seemingly independent of the number of the ben¬ 
zene rings. A careful analysis on the 2-ring case also 
showed that a possible superconducting phase with T c in 
the same range for K^CioHg. With these analysis, we 
are confident to conclude that the presence of the super¬ 
conducting phase with T c around 5 — 7 K is a common 
feature of aromatic hydrocarbons. In most cases, the 
low-T c phase was easily fabricated and more stable than 
the higher T c phase [7, 8, 10]. Considering the fact that 
the cation concentration has never been experimentally 
identified by using the reliable neutron scattering or sin¬ 
gle crystal X-ray diffraction techniques, all the reported 
doping levels are an ideal chemical estimation based on 
the synthesis procedure. Our phase stability examina¬ 
tion for potassium-doped benzene pins down the doping 
concentration of the exact two electrons for the stabile 
superconducting phase of these hydrocarbons with T c of 
5 - 7 K. 

Although a complete understanding of the mechanism 
of superconductivity in aromatic hydrocarbons has yet 
to achieve, the common feature of the 5 — 7 K super¬ 
conducting phase for all reported aromatic hydrocarbons 
may provide us some clues. Moreover, further investi¬ 
gations indicate that the electron-electron correlations 
seem to play an important role on superconductivity in 
these systems. This is quite similar to charge-transfer or¬ 
ganic superconductors and fullerides [5, 40-45] where the 
importance of electron correlations on superconductivity 
was emphasized. Magnetic measurements [9] on phenan- 
threne suggest the existence of local spin moments, and 
indicate similar interplay between magnetism and super¬ 
conductivity observed in fullerides [5, 40-45]. We indeed 
found that the electron-electron correlation becomes im¬ 


portant as the number of benzene rings increases. Multi¬ 
ring hydrocarbons are in favor of strong correlations and 
thus tend to have higher T c ’s. It is therefore believed 
that both the electron-phonon interaction and electron- 
electron correlations work together to enhance T c of the 
aromatic hydrocarbons with increasing the number of 
benzene rings. 

In summary, carrying out the first-principles calcula¬ 
tions, we have investigated the crystal structures and 
possible superconductivity in potassium doped solid ben¬ 
zene, KHg (x = 1,2,3). The K dopant results in 
the visible change of crystal configuration and drives the 
system from insulator to metal. K^CgHg is found to be 
stable due to its lowest formation energy as compared to 
other doping. Superconductivity in doped solid benzene 
with a T c of 6.2 K is predicted. This superconducting 
phase shares the common feature of all existing aromatic 
hydrocarbon superconductors. 
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